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Abstract 



We have searched for first generation scalar leptoquark (LQ) pairs in the 
e^+jets channel using pp collider data (/ Ldt « 115 pb -1 ) collected by the 
D0 experiment at the Fermilab Tevatron during 1992-96. The analysis yields 
no candidate events. We combine the results with those from the ee+jets 
and vv+]ets channels to obtain 95% confidence level (CL) upper limits on 
the LQ pair production cross section as a function of mass and of (3, the 
branching fraction to a charged lepton. Comparing with the next-todeading 
order theory, we set 95% CL lower limits on the LQ mass of 225, 204, and 
79 GeV/c 2 for (3 = 1, |, and 0, respectively. 



PACS numbers: 12.60.-i, 12.90.+b, 14.80.-j, 13.85.Rm 
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One of the remarkable features of the Standard Model (SM) is the symmetry between 
quarks and leptons that leads to cancellation of chiral anomalies and renders the SM renor- 
malizable. This symmetry might be explained by a more fundamental theory that relates 
quarks and leptons. Several extensions of the SM include leptoquarks (LQ), color-triplet 
bosons which carry both lepton (£) and quark (q) quantum numbers. The masses and 
coupling strengths of leptoquarks that couple to all three fermion generations are severely 
constrained by low energy experiments |2j and by HERA ||. Therefore, only LQ that cou- 
ple to a single generation can be light enough to be accessible at present accelerators. The 
excess of events at high Q 2 in e + p collisions reported || by the HI and ZEUS experiments 
at HERA, and its possible interpretation H as evidence for production of first generation 
scalar leptoquarks with a mass near 200 GeV/c 2 , have stimulated much interest in these 
particles. 

Leptoquarks would be dominantly pair-produced via strong interactions in pp collisions, 
independently of the unknown LQ-£-q Yukawa coupling. Each leptoquark would subse- 
quently decay into a lepton and a quark. For first generation leptoquarks, this leads to 
three possible final states: ee+jets, ez^+jets and vv+jets, with rates proportional to [3 2 , 
2(3(1-/3), and (1 — (3) 2 , respectively, where (3 denotes the branching fraction of a leptoquark 
to an electron and a quark (jet). The CDF and D0 |||| Collaborations have both 
searched for first generation scalar leptoquarks. The recent analyses ]9|,[7j of the ee+jets 
decay channel yielded 95% confidence level (CL) lower limits on the leptoquark mass of 
225 GeV/c 2 (D0) and 213 GeV/c 2 (CDF). In this Letter we present an analysis of the 
ez/+jets final state (which has maximum sensitivity at (3 — |), using 115 ±6 pb _1 of collider 
data collected at the Fermilab Tevatron at ^Js = 1.8 TeV during 1992-96. We also present 
a reinterpretation of our search [T(J for top squark pairs as a search for leptoquarks in the 



vv+jets decay channel. We combine the results from all three decay channels to obtain 
lower limits on the leptoquark mass as a function of (3. 

The D0 detector |j | consists of a central tracking system including a transition radiation 
detector, a uranium/liquid- argon calorimeter and a muon spectrometer. The data used in 
this analysis were collected with triggers which required the presence of an electromagnetic 
object, with or without jets and missing transverse energy ($ T ). The combined efficiency 
of the triggers is greater than 98% for LQ masses above 80 GeV/c 2 . Offline event selection 
requires: one electron with transverse energy Ej, > 20 GeV and pseudorapidity < 1.1 
or 1.5 < |?7| < 2.5; Ifi T > 30 GeV as the signature for a neutrino; and two or more jets 



reconstructed using a cone algorithm (cone radius 1Z = y (A0) 2 + (Ar?) 2 = 0.7, where is 

the azimuthal angle) with E° T > 20 GeV and \r]\ < 2.5. To reduce the effects of jet energy 
mismeasurements, we require the $ T vector to be separated from the jets by A0 > 0.25 
radians if $ T < 120 GeV. To suppress background from heavy quark events, we reject 
events that contain muons. These selection criteria define our base sample. 

An electron is identified by its pattern of energy deposition in the calorimeter, the pres- 
ence of a matching track in the central tracking detectors, and ionization in the central detec- 
tors. The efficiency for finding an electron is calculated to be (61±4)%, using Z(— > ee)+jets 
events which are similar in topology to the signal. All kinematic quantities in the event are 
calculated using the event vertex determined by the electron. 

We use the ISAJET pi Monte Carlo (MC) event generator to simulate LQ signal events 
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for masses (M LQ ) between 80 and 220 GeV/c 2 in 20 GeV/c 2 steps. The PYTHIA 0] MC 



program is used to simulate signal events with M LQ = 200 GeV/c 2 to study systematic 
errors arising from differences in modeling of the gluon radiation and parton fragmentation. 
The leptoquark production cross sections used are from recent next-to-leading order (NLO) 
calculations fl5| . The dominant VF+jets background is simulated using the VECBOS |T6|j 
event generator (with the HERWIG [1_7| program used for fragmenting the partons). The 
background from multijet events containing a jet which is misidentified as an electron, and 
with Ifjrp arising from the mismeasurement of jet energies, is modeled using multijet data. 
The probability for a jet to be misidentified as an electron (the fake probability) is estimated 
to be (3.50 ±0.35) x 10~ 4 . The background from tt decays into one or two electrons plus 
two or more jets, is simulated using the HERWIG MC program with a top quark mass of 
170 GeV/c 2 . All MC event samples are processed through the D0 detector simulation based 
on the GEANT |TJ|] package. 

In the base sample of 1094 events, we estimate the number of tt events to be 12 ± 4 



using the measured tt production cross section of 5.5 ± 1.8 pb [12]. The multijet background 
is estimated to be 75 ± 15 events, using a sample of events with three or more jets with 
]firp > 30 GeV. This is done by multiplying the fake probability by the number of ways 
the events satisfy the selection criteria with one of the jets passing the electron Ej- and 
f] requirements. After the estimated numbers of tt and multijet background events are 
subtracted, the number of events with transverse mass of the electron and neutrino (M T V ) 
below 110 GeV/c 2 is used to obtain an absolute normalization for the VF+jets background. 
This background is then largely eliminated by requiring M T U > 110 GeV/c 2 . After this 
cut, 14 events remain in the final data sample. The estimated background is 17.8 ±2.1 
events, of which 11.7 ± 1.8, 4.1±0.9, and 2.0±0.7 events are from VF±jets, multijets, and tt 
production, respectively. Leptoquark pair production would yield 24 events in this sample, 
if M L q = 120 GeV/c 2 and f3 = \. Assuming all 14 events to be signal, LQ production 
for masses below 120 GeV/c 2 can be excluded at the 95% CL for (3 = ~ with no further 
optimization. 

We have identified two additional variables that provide significant discrimination be- 
tween signal and the remaining background. They are the scalar transverse energy sum 
St = E T + Ej, + E T 2 + $ Tl where E T ia are the transverse energies of the two leading jets, 
and a mass variable ^(Mlq) = min( ^ MeJ | fL) ^ /LQ , ^^m^^ )■> where Mlq is an assumed LQ 
mass and M e j 1 2 are the invariant masses of the electron with the first and second leading 
jets. 

To find the optimal selection cuts, we adopt the criterion || of maximizing the MC signal 
efficiency for a fixed expected background of approximately 0.4 events. In the low mass 
range (Mlq < 120 GeV/c 2 ), where LQ production rates are high, requiring St > 400 GeV 
is sufficient. For Mlq > 120 GeV/c 2 , we use neural networks (NN) since they provide higher 
efficiency than an St cut alone. At each mass, Mlq, where we have generated MC events, we 
use a three layer feed- forward neural network [|B| with two inputs (St and ^(Mlq)), five 
hidden nodes, and one output (X ) nn(M L q)). We train each NN using simulated LQ events 
as the signal (with desired £> nn (M L q) = 1) and a mixture of VK±jets, multijet, and it events 
as background (with desired £> nn (M L q) = 0). Cuts on £>nn(M L q) that yield background 
estimates closest to the desired background are obtained by varying £>nn(M L q) in steps of 
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0.05. The background after the cut ranges between 0.29 ± 0.25 and 0.61 ± 0.27 events as 
shown in Table |. The errors in the background estimates include the effects of uncertainties 
in jet energy scale, fake probability and tt production cross section. The signal detection 
efficiencies calculated using simulated LQ events passing the selection requirements are also 
shown in Table |. The errors on the signal efficiencies include uncertainties in trigger and 
particle identification efficiencies, jet energy scale, effects of gluon radiation and parton 
fragmentation in the signal modeling, and finite MC statistics. No data events pass the 
cuts. 

To demonstrate that the backgrounds are reliably modeled, comparisons of the data and 
combined background in the variables M T U and Pnn(180) are shown in Fig. [I] for the base 
sample. 

Figures |2| (a)-(c) show the 2-dimensional distributions of ^(180) vs. St for simulated 
LQ signal events with Mlq = 180 GeV/c 2 , the combined background, and data. The 
contours corresponding to constant values of £>nn(180) demonstrate the separation achieved 
between signal and background. The distribution of £>nn(180) for data is compared with 
the predicted distributions for background and signal in Fig. ^| (d). It is clear that the data 
are described well by background alone. The highest £>nn(180) observed in the final data 
sample is 0.79. 

Using Bayesian statistics, we obtain a 95% CL upper limit on the leptoquark pair pro- 
duction cross section for (3 = | as a function of leptoquark mass. The results are shown in 
Table |[ The statistical and systematic uncertainties in the efficiency, the integrated lumi- 
nosity, and the background estimation are included in the limit calculation with Gaussian 
prior probabilities. The measured 95% CL cross section upper limits for LQ pair production, 
corrected for the branching ratio with [3 = |, for various LQ masses are plotted in Fig. |3| 
together with the NLO calculations The intersection of the limit curve with the lower 
edge of the theory band (renormalization scale \i = 2Mlq) is at 0.19 pb, leading to a 95% 
CL lower limit on the LQ mass of 175 GeV/c 2 . 

An analysis of the vv + jets channel is accomplished by making use of our published 



search (with J Ldt « 7.4 pb ) for the supersymmetric partner of the top quark [ T0| . Three 



events survive the selection criteria ($ T > 40 GeV, 2 jets with E T > 30 GeV, and no isolated 
electrons or muons) consistent with the estimated background of 3.5 ±1.2 events, mainly 
from iy/Z+jets production. The efficiencies of the event selection for Mlq= 60, 80, and 
100 GeV/c 2 are calculated to be 1.1%, 2.2%, and 3.9%, respectively, using signal MC events 
generated with the ISAJET generator and processed through the detector simulation based 



on GEANT. The systematic errors in the signal acceptance are calculated as in Ref. |T0 
This analysis yields the limit M LQ >79 GeV/c 2 at the 95% CL for =0. 

Combining the ee+jets, ei^+jets, and vv+jets channels, we calculate 95% CL upper limits 
on the LQ pair production cross section as a function of LQ mass for various values of f3. 
These cross section limits for j3 = | (shown in Fig. [$]), when compared with NLO theory, 
yield a 95% CL lower limit on the LQ mass of 204 GeV/c 2 . The lower limits on the LQ mass 
derived as a function of (3, from all three channels combined, as well as from the individual 
channels, are shown in Fig. f|. These results can also be used to set limits on pair production 
of any heavy scalar particle decaying into a lepton and a quark, in a variety of models. 

In conclusion, we have presented a search for first generation scalar leptoquark pairs in 
the e^+jets decay channel. Combining the results with those from the ee+jets and z/^+jets 
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channels, we exclude leptoquarks with mass below 225 GeV/c 2 for (3 = 1, 204 GeV/c 2 for 
= 5, and 79 GeV/c 2 for /3 = 0, at the 95% CL. Our results exclude (at the 95% CL) the 
interpretation of the HERA high Q 2 event excess via s-channel scalar LQ production with 
LQ mass below 200 GeV/c 2 for values of j3 > 0.4 and significantly restrict new LQ models 



containing additional fermions [21 
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FIG. 1. (a) M?p u and (b) 2?nn(180) distributions for data (points) and background (histograms), 
after all cuts except on Mjf and X>nn(180). The arrow in plot (a) shows a cut on Mjf, as described 
in the text. 
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FIG. 2. Distributions of ^-(180) vs. St for (a) predicted background, (b) simulated LQ 
events (Mlq = 180 GeV/c 2 ), and (c) data, after all cuts except that on X>nn(180). The contours 
correspond to Pnn(180) = 0.75, 0.85, and 0.95. The box area is proportional to the number 
of events in the bin, with the total number of events normalized to 115 pb" 1 . Plot (d) shows 
distributions of X>nn(180) for data (solid circles), background (open histogram) and expected LQ 
signal for Mlq = 180 GeV/c 2 (hatched histogram). The arrow in plot (d) shows the chosen cut 
on X>nn(180), as described in the text. 
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TABLE I. Signal detection efficiencies, estimated backgrounds and measured 95% CL upper 
limits on the production cross section from the ez^+jets channel analysis. The NLO cross sections 
(with /x = 2Mlq) from Ref. [15] times 2/3(1 — (3) = ^ for (3 = ~ are also shown. 



Leptoquark 


Signal 


Estimated 


95% CL 


NLO 


Mass 


Efficiency 


Background 


Upper Limit 


Theory 


(GeV/c 2 ) 


(%) 


(Events) 


(pb) 


(Pb) 


80 


0.3 ±0.1 


0.60 ±0.27 


10.88 


17.98 


100 


1.2 ±0.2 


0.60 ±0.27 


2.59 


5.34 


120 


2.5 ±0.3 


0.60 ±0.27 


1.15 


1.90 


140 


6.7 ± 1.0 


0.54 ±0.25 


0.43 


0.77 


160 


10.9 ±1.2 


0.61 ±0.27 


0.25 


0.34 


180 


14.7 ± 1.2 


0.29 ±0.25 


0.18 


0.16 


200 


19.4 ±1.7 


0.43 ± 0.27 


0.14 


0.08 


220 


21.5 ± 1.7 
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FIG. 3. Measured 95% CL upper limits on the leptoquark pair production cross section (see 
text) in the eu + jets channel (circles) and all three channels combined (triangles) for /3 = |. Also 
shown are the NLO calculations of Ref. [15| where the central line corresponds to /x = Mlq, and 
the lower and upper lines to (i = 2Mlq and \i = ^Mlq, respectively. 
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FIG. 4. Lower limits on the first generation scalar leptoquark mass as a function of /3, based 
on searches in all three possible decay channels for LQ pairs. Limits from LEP |2C] and from our 
previous analysis M of 1992-93 data are also shown. The shaded area is excluded at 95% CL. 
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